A simple model for DNA denaturation 
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Abstract 

Following Poland and Scheraga, we consider a simplified model for the de- 
naturation transition of DNA. The two strands are modeled as interacting 
polymer chains. The attractive interactions, which mimic the pairing be- 
tween the four bases, are reduced to a single short range binding term. Fur- 
thermore, base-pair misalignments are forbidden, implying that this binding 
term exists only for corresponding (same curvilinear abscissae) monomers of 
the two chains. We take into account the excluded volume repulsion between 
monomers of the two chains, but neglect intra-chain repulsion. We find that 
the excluded volume term generates an effective repulsive interaction between 
the chains, which decays as l/r d ~ 2 . Due to this long-range repulsion between 
the chains, the denaturation transition is first order in any dimension, in 
agreement with previous studies. 

Submitted for publication to: "Europhys. Lett." Saclay, SPhT/01-003 



1 



The problem of the denaturation of DNA has been of physical interest for a long time 
]1].|2|], and can be roughly characterized as the adsorption (or more accurately, desorption) 
transition of one strand onto (from) the other. Experiments indeed show that the fraction 
of bound base pairs exhibits, as a function of temperature, a series of sharp jumps, linked 
to the local melting of these heterogeneous paired structures. This is in contrast with many 
theories |],|2|], which yield a continuous denaturation transition below four dimensions. 

The full problem is rather complicated since different binding energies are associated 
with different bases contents (such as different AT and GC frequencies). 

Moreover, since the denaturation transition decreases the stiffness of a DNA molecule, 
the backbone elasticity may play an important role in the transition j3|. Along these lines, 
the helical structure of the molecule has been also considered [0 

To disentangle polymeric from disorder effects, simplifying assumptions have to be made. 
In this work, we follow a time honored approach where (i) disorder is neglected (homoge- 
neous DNA approximation) (ii) base-pair mismatch is forbidden (pairing is possible only 
between corresponding monomers along the two chains), (hi) we consider excluded volume 
interactions only between the two chains. 

Experimentally, the chains are not very long (N « 1000) and have a fairly large persis- 
tence length (l p ~ 20 pairs). This suggests that for typical distances of the chains smaller 
than several persistence lengths, excluded volume effects between monomers in the same 
chain will be negligible compared to excluded volume effects between monomers of the two 
chains. 

Closely related work can be found in references ||,^| where the existence of a first order 
transition in three dimensions has been clearly linked, via numerical simulations or scaling 
arguments, to the short range (excluded volume) repulsion between the chains (by first 
order, we mean a phase transition with a finite latent heat and/or some finite length scales). 
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In our model, the partition function of the two chains of length N reads : 



xexp^-gj ds J ds' 5{r x {s) — f^s'))! (1) 

where d is the space dimension, a is the Kuhn length of the monomers, j3 is the inverse 
temperature, g is the excluded volume parameter and v(f\ — r 2 ) is the short range binding 
potential for monomers s of chain 1 and 2. In the following, we will model this interaction 
by an attractive (d dimensional) spherical well potential, of radius r and depth — Vq. 

An approximate treatment of the excluded volume interaction will be given first. Then, 
we study the unbinding phase transition in the framework of a quantum analogy @. We 
finally compare our results with references |||6| . 

It is clear from equation ([[]), that one should focus on the relative coordinate (ri(s) — 
r*2(s)) of the two chain system to study a possible adsorption transition. Introducing the 
new coordinates 

r(s) = m±m (2) 

f(s) = fi(s) -f 2 (s) (3) 
and performing the change of variable in ([]]), we have 

Z = |pf( S )exp(-A j" & -pj"dsv(r(s))-W({r(s)})\ (4) 



where 

-W({r(s)}) 



d f N ^ (dR\ 2 „ r»j_ Fj.J*,^ s,^ , r(s)+f(s'y 



JvR(s)exp \—^f ds ^ -gj ds j ds'5 (r(s) - R(s') + 



2 

(5) 



In the denatured phase, the two chains are far apart and interact very weakly; it is 
therefore justified to treat this interaction in a perturbative manner. Expanding the effective 
potential W to first order in g, we get 



-W({r(s)}) 



exp l-g J q N ds J q N ds' (5 (r(s) - R(s') + ^(f(s) + *V)))) ) (6) 



where the brackets stand for 

7 _ AT / in\ 



(A) =i/M W Xe X p^f ds [f) j < 7 > 

and V is the volume of the system. 
We obtain 

which is valid for dimensions d > 2. This expression is asymptotically exact in the denatured 
phase, where very few monomers of the two chains come close together. 
Setting S = (s + s')/2 and a = s — s' in equation <^), we get for large N 

N jo f + ™ j ( d ( d{r{S + a/2) + r{S -a/2)) 2 ^ 







wm.)})*9l_ iSf ^\-—A exp ^ J (») 



A simple scaling argument shows that this effective interaction must vanish at large 
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relative separation. The range of |<r| which contributes to the integral is given by |c| 
d(f(S + cr/2) + f(S — a/2)) 2 /4a 2 . Thus, for small separations where the interaction is 
sizeable, only small a contribute significantly to the integral and we may expand f(S + 
cr/2) + f(S — cr/2) in powers of a, and obtain 

v r+00 ( d \ d/2 ( 



W({r(s)})~g dS exp 

JO J -00 V 7T <T d I 



a 2 \o\ 



N a d 



, ds JT^ (10) 

\r(s)\ 



where 



*< = M-) n- 2 -D- 2 (id 



We see that the effect of the inter-chain excluded volume interaction is to generate a 
"Coulomb" -like repulsive term in the relative coordinate. This entropic repulsion is familiar 



in other interacting lines problems PJTOfl. Its long range nature is due to the large lateral 
excursions of the chains. 

An alternative derivation of this term can be obtained from equation ([I]). The use of the 
identity 

5(n(s) - f 2 (s')) = I C *<*«-*C» (12) 
allows us to rewrite the excluded volume term of equation (0) as 

- j N ds j N ds' f -^-( e *((ft(«HM)+(ftM-ft(0)) + e ^-((^W-ft(«)+(ftW)-ft(«')))) (13) 
2 Jo Jo J (2ir) d 

where we have explicitly introduced, in a symmetric way, the relative coordinate of the 
chains (rl(s) — r 2 (s)). 



Each term of equation (|T^) is a product over single chain and relative coordinate confor- 
mations. For each contact of the chains, it seems natural to average the single chain term. 
Defining this average by G(k, s, s') = ^e tk ^ r ^~ r ^ s '^), the excluded volume term reads 

>E"£«I^ a $<-<<>**" Hm (14) 

Since we have neglected the intrachain excluded volume, we have 

G(k,s,s') = e-^ s - s ' 1 (15) 



Putting everything together, we get back equations (10,11). 



We thus can write the partition function (Q) for the relative coordinate of the two chains 

as 

For long chains, we may now use a well known quantum mechanical analogy |7j and 
study the ground state of the Hamiltonian 



This ground state (wavefunction *5/ , energy E ) is given by the Schrodinger equation 

H* = E * (18) 

The free energy per monomer then reads F = E Q T. 

We now discuss this equation as a function of space dimensionality. 

• 2 < d < 4 

We discuss the typical case d = 3. The potential is radial, and thus the ground state 
will be spherically symmetric. 

In radial coordinates, equation (|i"8|) becomes 



2 

~ W u o +Pv(r)u + — u = E u (19) 

O T 

where Uq{t) = r\l>o( r ) and «3 is given by (|TTJ). 

The unbinding transition [^,3 occurs when T = T c and E = 0. As the binding 
potential vanishes for r > tq, we have in this region 

-^< + -«o = (20) 

The solution of this equation which vanishes at infinity is given by a Bessel function 

Mr) = f x K 1 (^) (21) 

with 

u (r) ~ ^— ^ exp ( K ~)J~~^J for r ^ oo (22) 
The (finite) correlation length A is given by 

implying that the transition is first-order. This is due to the strongly repulsive barrier 
generated by the excluded- volume. The transition temperature T c is obtained by 
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matching the above solution to the solution inside the attractive region (v(r) = — Vq 
for r < r ). The calculations are rather tedious and will not be given here. It is 
however interesting to derive an estimate of the latent heat at the transition. We 



therefore rewrite equation (|T9|) as 



2 

\u n Q + f3 c v{r) u +(p- f3 c )v{r) u + —u = E u (24) 



Close to T c , we do first order perturbation theory and get 

E = (u 0c (r)\(p-(3 c )v(r)\uoc(r)) (25) 

where uo c (r) denotes the full solution of equation fll9|) at the critical temperature. 
Since F ~ T C E , we finally obtain for the latent heat L (L = T c {^)t c ) 

L = -(u 0c (r)\v(r)\u 0c (r)} (26) 



A simpler procedure [0] replaces the attractive region by an attractive 5— potential of 
strength 7 at r = r leading to 

F = -7 i^P) *Z(r ) for T - T~ (27) 
= for T > T c (28) 

and a latent heat L = juKro), where uo(r ) is given in equation fl21~|) . 
• d > 4 

The ground state \l/o satisfies the radial equation 

a 2 . d 2 d — 1 d w / _ , . «d\ T 

The unbinding point occurs for E = . It is easily seen that at large distances, the 
"Coulomb" potential is negligible, and thus the wave-function decays as a power-law 
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* (r) ~ for r -> oo (30) 

This wave-function is normalizable for d > 4, and thus, the transition is again first- 
order. This is in agreement with the fact that in d > 4, the excluded-volume effect 
should be irrelevant and thus one should recover the results of the Gaussian adsorption 
problem. According to the Poland- Scheraga theory [|l|, the denaturation transition is 
first-order if the exponent 9 describing the probability of first return to the origin of 
the chain is larger than 2. For a Gaussian chain in d > 4, this is indeed the case since 

e = i 

However, this wave-function is much less localized than below four dimensions: only 
moments of r of order strictly smaller than (d — 4) are finite. 

d = 4 

This case is marginal, since the Laplacian operator and the " Coulomb" potential have 
the same scaling dimension. At large distance, the wave-function satisfies 

fl2 ( ^ 2 + ^ ^ )\& + a4 \& (31) 
4 dr 2 r dr r 2 



The wave-function behaves as a power-law 



*o ~ p (32) 



with 



v 



+ (33) 



The wave-function is normalizable, but again with an infinite correlation length and only 
a finite number of moments of r are defined. 

We have studied a simple model for the denaturation of DNA, where we have given 
an asymptotically exact treatment of the interchain excluded volume interaction. We find 
that the first order character of the denaturation transition stems from the existence of 
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this long range interaction, in broad agreement with recent work 0,^]. The difference in 
our results might be due to our neglect of the intrachain excluded volume. Within their 
scaling approximation, the latter authors find that the first order transition is associated 
with a large loop probability distribution, whereas our approximation yields an exponentially 
bound state for d = 3. Future work is needed to clarify this point. Sequence heterogeneity 
11 1, and chain stiffness |l^JT3[| raise further questions. 
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